Pentacene multilayers on Ag(111) surface 
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Abstract 

The structural profiles and electronic properties of pentacene (C22H14) multilayers on 
Ag(l 11) surface has been studied within the density functional theory (DFT) framework. We 
have performed first-principle total energy calculations based on the projector augmented wave 
(PAW) method to investigate the initial growth patterns of pentacene (Pn) on Ag(l 11) surface. 
In its bulk phase, pentacene crystallizes with a triclinic symmetry while a thin film phase hav- 
ing an orthorhombic unit cell is energetically less favorable by 0.12 eV/cell. Pentacene prefers 
to stay planar on Ag(l 11) surface and aUgns perfectly along silver rows without any molecular 
deformation at a height of 3.9 A. At one monolayer (ML) coverage the separation between 
the molecular layer and the surface plane extends to 4.1 A due to intermolecular interactions 
weakening surface-pentacene attraction. While the first ML remains flat, the molecules on a 
second full pentacene layer deposited on the surface rearrange so that they become skewed 
with respect to each other. This adsorption mode is energetically more preferable than the one 
for which the molecules form a flat pentacene layer by an energy difference similar to that 
obtained for bulk and thin film phases. Moreover, as new layers added, pentacenes assemble 
to maintain this tilting for 3 and 4 ML similar to its bulk phase while the contact layer always 
remains planar. Therefore, our calculations indicate bulk-like initial stages for the growth pat- 
tern. 
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Introduction 



Due to its use in thin film transistor (TFT) applications pentacene is continuing to enjoy being 
the subject of extensive research. Since most TFT's employ Si02 as the dielectric layer, initial 
studies on pentacene focused mainly on gaining a thorough understanding of structural ^"^ and 
electronic^"- properties of thin films of this molecule on Si02 surfaces and in turn achieving the 
best device performance by optimizing-"— these properties. As a result of this heavy research 
effort in the last decade, fabricating pentacene TFT's with hole mobilities of more than 1 cm^/Vs 
have become an almost routine process.— Nevertheless there are still, fundamental issues to be 
resolved such as the dependence of the charge mobility on the film thickness^ and areas open to 
improvement like modification of the substrate surfaces with buffer layers- or the pentacene itself 
with functional groups.-*^ 

Another very critical but relatively less well understood subject is the growth mechanism of 
pentacene films on metal substrates. Understanding the pentacene film growth on gold and silver 
surfaces is particularly important since these metals constitute the electrode material in most TFT's 
and the device performance is directly related to the charge transfer efficiency between the elec- 
trodes and the organic film. Though both experimental—"— and theoretical-^^'^i^'^ research in this 
field has been recently intensified, there are contradictory results in the literature and the growth 
modes of pentacene thin films on Au(l 1 1) and Ag(l 11) surfaces are continuing to be a matter of 
debate. This is mostly due to the relatively strong (when compared with Si02) interaction of the 
metal surfaces with the pentacene molecule. As a result of this strong interaction pentacene adopts 
many different monolayer and multilayer phases on metal surfaces which are energetically and 
structurally very close to each other. For example on Au(l 1 1) several different low density mono- 
layer phases and an identical full coverage phase have been reported by different groups.— "—>^ 
However, in case of the multilayers while Kang et al.— report a layer by layer growth of ly- 
ing down pentacene molecules, Beemink et al. report strong dewetting starting from the second 
layer and growth of bulk like pentacene crystals. For pentacene films on Ag(l 11) surfaces, while 
Eremtchenko et al.— and Dougherty et al.^ report a bilayer film formation, where an ordered 
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(second) layer (which follows the symmetry of the Ag(l 11) surface) forms on top of a disordered 
(2D gas phase) first layer at room temperature, Kafer et al.^^ report the formation of bulk-like 
pentacene structures immediately after the first monolayer. So on both surfaces the growth mech- 
anism of pentacene films is still not completely clear. In addition, if the above mentioned 2D gas 
phase mechanism for Ag(lll) is really correct then questions like "Why does pentacene behave 
completely differently on seemingly similar surfaces, Ag(l 1 1) and Au(l 1 1)?" and "How does the 
symmetry of the substrate affect the bilayer film structure?" arise. 

In spite of this richness of experimental studies and points in need of clarification, pentacene 
films on Ag(lll) or Au(lll) surfaces have not, yet, been studied theoretically. Theoretical work 
regarding pentacene films were mostly performed on other metal surfaces, such as Cu(OOl),— 
Cu(llO),^ Ag(110),-2^ at semi-empiric level and Al(100),35 Cu(100),3^ Cu(119),-22 Fe(lOO),^ 
Au(OOl)-^^'-^^ at DFT level. These studies in general addressed two points concerning the first 
layer of pentacene film/molecule: (1) Determination of the most stable adsorption site/geometry, 
and (2) determination of the strength of electronic interaction/coupling between the substrate and 
the molecule. In these studies either the most stable configuration was found to be pentacene ly- 
ing flat on the surface or the calculations were started with this assumption. In terms of the 
electronic interactions, the DFT studies performed using GGA functionals found considerable 
aromatic-;r-system metal substrate interaction—"^ on Cu surfaces, hinting at chemi sorption. On 
Au(lll)— and Al(lOO)— however, while LDA functionals resulted in strong interactions, in the 
form of broadening and splitting of ;r-molecular orbitals, GGA functionals are reported to result 
in much weaker interactions, in accord with a physisorption mechanism. Theoretical studies con- 
cerning the further stages of pentacene film growth on metal surfaces, however, like second layer 
structure/energetics or the thin film crystal/electronic structure, are very few and at semi-empiric 
level. Instead, theoretical works regarding pentacene films are primarily focused on the elec- 
tronic structure of different pentacene polymorphs observed mainly on SiOi surfaces, one being 
the famous "thin film phase" .■^^^'^ 

Hence, a theoretical study of growth mechanism and electronic properties of pentacene films 
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on Ag(lll) and Au(lll) may, (1) help resolve the experimental contradictions mentioned above 
and (2) fill a gap in the theoretical literature and enable a comparison of these systems with other 
pentacene films. As a first attempt to this end, here we present the results of our work on the 
structural and electronic properties of monolayer and multilayer films of pentacene on Ag(l 11) at 
DFT level. First we discuss full coverage monolayer film in the light of experimental results. We 
compare the adsorption geometries and the corresponding density of states we found with the ex- 
perimental results reported so far. Then we present our results regarding two and three monolayers 
of pentacene film and discuss how the crystal and electronic structure of the Ag(lll)/pentacene 
interface and the film evolves with coverage. We conclude with an overall summary and discussion 
of the results. 

Method 

We performed total energy density functional theory (DFT) calculations using the projector- augmented 
wave (PAW) method— within the generalized gradient approximation (GGA) by employing the 
Perdew-Burke-Ernzerhoff (PBE)— exchange-correlation energy functional as implemented in the 
Vienna Ab Initio Simulation Package (VASP).— 

For consistency, we used a kinetic energy cutoff of 370 eV for the plane wave expansion of 
single particle wavefunctions in all the calculations. Electronic ground states has been determined 
by requiring a total-energy convergence up to a tolerance value smaller than 0.1 meV. We used a 
conjugate-gradient algorithm, in all geometry optimization calculations, based on the reduction in 
the Hellman-Feynman forces on each constituent atom to less than 10 meV/A. 

We examined two previously known polymorphs of Pn lattice: the bulk— and the thin film- 
phases. Bulk phase corresponds to a triclinic unit cell which contains two €22^14 formulae with a 
set of parameters: a = 7.90 A, ^ = 6.06 A, c= 16.01 A, a = 101.9°, /3 = 112.6°, and 7= 85.8°.-^ 
Our calculated values of a = 7.90 A, b = 6.06 A, c = 16.01 A, a = 102.0°, /3 = 112.6°, and 
Y— 85.5° shows a very good agrement with the experimental results of Campbell et al. For the 
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thin film phase, we obtained an orthorhombic unit cell with a = l A2 A,b = 5.87 A, and c = 16.21 
A. These parameters agree with Parisse et aUs^ theoretical results of a = 1 .60 A, b = 5.90 A, 
and c = 15.43 A except for the last one, corresponding to the longitudinal size of the unit cell. 
We calculated the length of an isolated molecule to be 14.5 A. Therefore, the disagreement, in this 
phase, stems from the Pn-Pn separation in the multilayers. 

In order to compare the predictions of different DFT functional with the experimental results 
for key structural and electronic parameters like the lattice constants and binding energies, we re- 
peated our calculations using PW91— parametrization within both generalized gradient and local 
density approximations (LDA). We calculated the lattice parameter of the silver ccp bulk struc- 
ture in Fm3m symmetry group to be 4.00, 4.15, and 4.17 A using LDA-PW91, GGA-PW91, and 
GGA-PBE functionals respectively. These compare well with the experimental^^ value of 4.09 
A, slightly better than the previous theoretical^*^ result of 4.20 A. The Ag(lll) surface has been 
modeled in a four-layer slab geometry separated from their periodic images by ~ 15 A of a vacuum 
space and 3x6x1 grid for the cases of IML to 4ML deposition, whereas in the cases of isolated Pn 
adsorption a larger cell is needed and therefore a three-layer- slab geometry for the Ag(l 1 1) surface 
and 1x2x1 grid was used. For this metallic system, the number of layers has been found to be 
sufficiently large to represent the Ag(lll) surface structure such that the geometry optimization 
calculations do not disturb the subsurface layer atoms from their bulk lattice positions. 

Results and Discussion 

Single isolated pentacene on Ag(lll) 

In order to study the formation of ordered Pn layers on Ag(lll) we first considered a single Pn 
on the surface. Isolation of the molecules was achieved by using a 8x5 silver surface unit cell 
which sets 9.5 A tip-to-tip and 7.9 A side-to-side separations between the periodic images. We 
determined the minimum energy Pn/Ag(lll) geometry by investigating all possible adsorption 
sites with a number of orientations, compatible with the lattice symmetry, at each site as shown 
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Figure 1: Single isolated pentacene on different adsorption sites of Ag(lll) surface. Planar Pn 
adsorption with the central carbon ring on top of a silver atom aligning parallel to one of the lattice 
directions is abbreviated as "Top-0" (a). Top-30 in (b) refers to the adsorption at the top site with 
Pn major axis making an angle of 30° with any of the silver rows. HoUow-O (c) and Hollow-30 (d) 
follow the same molecular alignments as the top cases, but centered on a triangle whose corners 
defined by Ag atoms, i.e. at the hollow site. Bridge-0 (e), Bridge-30 (f), Bridge-60 (g), and Bridge- 
90 (h), describe the cases where central ring of Pn lies over an Ag-Ag bond making the referred 
angles with any of the lattice lines. The minimum energy geometry, Bridge-60 (g), is depicted in 
ball-and-stick fashion while the others are all shown in sticks only. 



7 



in Figure 1 (Labeling conventions are described in figure caption.) In addition to these planar 
cases we also investigated the possibility of standing-up adsorption configurations which appeared 
to be around 0.15 eV less favorable than the planar ones. In geometry optimization calculations 
Pn develops a weak interaction with Ag(lll) wherever it is initially placed on the surface. In 
fact, as presented in ITable 1[ the comparison of the total energies of these adsorption cases show 
differences which are no greater than 36 meV from each other. The flatness of the potential energy 
surface (PES) is indicated by the existence of such small barriers which might make Pn diffusion 
over the surface possible in agreement with the experimental observations that the contact layer Pn 
molecules are mobile at the Pn/Ag(lll) interface.—"^ Similarly, during image acquisition STM 
tip has been observed to drag Pn molecules which are physisorbed on Au(l 1 1).— 
Table 1: Calculated values for geometrical and electronic structure of Pn/Ag(lll) systems shown 
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Figure 1 The lateral height of isolated Pn molecule from the Ag(lll) surface in A, the 



binding energy E}, and the relative total energy Ej in eV. 





dz 




Ej 


Top-0 


3.90 


-0.125 


0.030 


Top-30 


3.89 


-0.119 


0.036 


HoUow-O 


3.88 


-0.147 


0.008 


Hollow-30 


3.87 


-0.128 


0.027 


Bridge-0 


3.87 


-0.124 


0.031 


Bridge-30 


3.88 


-0.124 


0.031 


Bridge-60 


3.87 


-0.155 


0.000 


Bridge-90 


3.88 


-0.129 


0.026 



The adsorption configurations (in Figure 1| ) where the isolated pentacene follows the lattice 



symmetry so that the molecular charge density matches better with the surface charge density of 
silver rows are energetically more preferable. As a result, the total energy of the bridge-60 is 
smaller from that of the hoUow-O by only 8 meV. This also indicates that the flatness of the PES is 
relatively more pronounced along the lattice directions. 

Single isolated pentacene molecule finds its minimum energy configuration at the bridge-60 



position as depicted in ball-and-stick form in Figure 1 1. For this adsorption geometry, it is almost 



flat with a negligible bending at a height of 3.87 A which gives a weak binding energy of —0.155 
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eV. In geometry optimization calculations, for all possible initial adsorption configurations, both 
GGA functional predict a weak interaction between the Pn molecule and the Ag(lll) surface 
where LDA overbinds. (see [Tabled . In particular, GGA-PBE predicts that an isolated pentacene 
with a tilt about 15 degrees off the surface plane is only 3 meV unfavorable than the lowest energy 
flat geometry. This barrier is so small that the tilted pentacene does not relax back to planar 
geometry. 

Although GGA functionals result in a weak pentacene-silver interaction, the degree of this 
weakness is overestimated. Since, pure DFT results depend on the choice of the exchange- 
correlation functional, a hybrid-DFT with corrected exchange with dispersive interaction energy 
would give an improved description of the binding characteristics of such a weakly bound sys- 
tem.— 

Table 2: Calculated values for electronic and geometrical structure of Ag and Pn/Ag(l 11) systems 
in different exchange-correlation functionals. Lattice parameter of Ag(l 11) slab flAg in A, lateral 
heights (A) and the binding energies £/, (eV) of isolated and 1 ML Pn on the Ag(l 1 1) surface. 

clean slab isolated Pn 1 ML Pn 

fiiAg d.^ Kjj El-, 

LDA-PW91 4.000 2.46 -1.925 2.48 -1.753 
GGA-PW91 4.145 3.69 -0.234 3.94 -0.093 
GGA-PBE 4.174 3.87 -0.155 4.12 -0.078 



Full monolayer 

Full monolayer has been derived from the previously optimized bridge-60 configuration using an 
experimentally observed 6x3 silver surface unit cell.— "^2*2^ The relaxed geometry of this contact 



layer is shown in Figure 3 i-c. The molecules follow the surface symmetry and are aligned parallel 
to the silver rows. The distance between Pn and Ag(lll) at the interface varies with varying 
intermolecular interaction strengths. For instance, in the case of GGA-PBE, an isolated Pn stays 
3.9 A above the surface while this value extends to 4. 1 A in the case of 1 ML coverage as presented 
in lTable 21 Corresponding binding energy at IML is calculated to be 0.08 eV for GGA-PBE. This 
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shows a weak binding similar to the experimental observations.— '"■2^*2^ GGA-PW9 1 gives a slightly 
better description of both the Ag substrate and Pn contact layer geometries relative to GGA-PBE 
results. We also calculated the supercell total energies as a function of Pn-Ag(lll) distance as 



shown in Figure 2 for isolated and full monolayer cases. LDA incorrectly gives strong binding 
for this system since the charge densities are well localized around atoms leaving nearly empty 
interatomic regions. 
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Figure 2: Binding energy versus pentacene height on Ag(l 11) relative to the minimum energy po- 
sition (bridge-60), calculated with different exchange-correlation functionals and approximations 
both for a single isolated Pn and for IML Pn coverage. 



Our DFT calculations show that a small tilt angle of the contact layer will not yield a significant 
increase in the total energy with respect to flat geometry. This is due to the overestimated weakness 
of pentacene-silver interaction. Based on these DFT results we can not reject the possibility of an 
average tilt at the IML as well as the isolated single pentacene case depending on the experimental 
conditions. However, our calculations do not suggest a strong binding between the pentacene layer 
and the surface. Therefore, our calculations indicate a flat 1 ML physisorption rather than a tilted 
chemisorbed Pn layer which was concluded by Kafer et al. based on their NEXAFS and thermal 
desorption signatures. 
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Figure 3: Side views (a, b) and top views (c, d) of (IML, 4ML) pentacene, respectively, adsorbed 
on the Ag(lll) surface. No perspective depth is used in the production of the figures, hence the 
actual tilts are seen in (b) and (d) for the 4ML case. 

Multilayers 

We considered two different initial geometries for the second pentacene layer on the already opti- 
mized IML Pn/Ag(l 1 1) structure. For the first case, the second layer molecules are flat on the first 
monolayer where molecular axes follow the surface symmetry. In addition, a second layer pen- 
tacene stays above in between the two molecules underneath. The second case initial structure is 
the same as the first one except the molecular planes of only the second layer pentacenes are tilted 
around their long axes as observed in experimental studies.— "-^^^Geometry optimization calcula- 
tions resulted in a very small difference of 0.046 eV in the total energies in favor of the latter case 
in which second layer molecules slightly misaligned from the surface lattice direction by 6.4° in 
addition to the molecular plane tilting of 18° as presented in lTable 3l The smallness of the energy 
difference between the two cases can be addressed to the energy difference between the different 
phases of pentacene. For instance, we calculated the difference in the total energies between the 
bulk and the thin film phases of pentacene to be 0.12 eV for a cell having two molecular formulae 
units. 
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Figure 4: The charge density contour plot of 4 ML pentacene-Ag(l 11) interface on a plane normal 
to the surface and cut through a silver row. 

For the 2ML structure, the tilted layer (the second one) relaxes to a height of 3.6 A above the 
contact layer which is separated from the silver surface by 3.8 A. Resulting height of the first tilted 
pentacene layer from silver surface becomes 7.4 A. In the case of 3ML and 4ML structures this 

o o 

height converges to 7.2 A which is slightly lower than the experimental value of 7.8 A reported by 
Danisman et al. . 

Table 3: Calculated values for geometrical structure of Pn/Ag(l 1 1) systems where is the distance 
between first-layer pentacene and Ag(lll) surface, dn-m is the distance between nth and mth Pn 
layers (all in A). 9 and a„ are the tilt angles of the nth layer molecules about the (lll)-axis and 
about their major axes, respectively. 

d^ di-2 d2-3 di,-4 6 a2 oci, 

isolated 3.9 - - - 0.0 - - - 

IML 4.1 - - - 0.0 - - - 

2ML 3.8 3.6 - - 6.4 18 - - 

3ML 3.7 3.5 3.6 - 6.0 25 -17 - 

4ML 3.7 3.5 3.5 3.6 4.0 20 -23 15 



The flat and tilted pentacene configurations above the first layer has also been considered for 
the 3ML and 4ML initial structures. The difference in the total energies has been calculated to be 
0.127 eV and 0.125 eV in favor of the tilted molecules on the first layer for 3ML and 4ML cases, 
respectively. Therefore, bulk- like pentacene formation on Ag(l 1 1) surface above the contact layer 
is more preferable than flat lying multilayers. 
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In the case of multilayers, the separation of the top layer from the layer underneath is 3.6 A and 
all the inner interlayer distances become 3.5 A while the height of the contact layer converges to a 
value of 3.7 A after the third ML. In addition, our calculations show a decrease in the misalignment 
of pentacenes above the contact layer as new layers deposited on Ag(lll) surface up to 4ML. 
All multilayer geometries can be seen through Figure 3p -d, which only show views of the 4ML 
structure, since the interlayer distances and angles do not change significantly as new layers added. 
We also present the corresponding electronic density contour plot of pentacene-Ag(l 1 1) interface 



at 4ML in Figure 4 which indicates charge localization around pentacene molecules and in the 
substrate suggesting weak binding of the contact layer. 

Consequently, the thin film formation starting from the second ML indicates that the inter- 
molecular interactions are relatively stronger than the pentacene-silver interaction. Substantia- 
tively, the misalignment of pentacenes above the contact layer can also be attributed to the weak- 
ness of Pn-Ag(lll) interaction. In fact, this observation is in parallel with the experimental ob- 
servations of Danisman et al.^-^ whereupon desorption of the pentacene multilayers on a stepped 
Ag(lll) surface a new monolayer phase was observed. In addition, tilted second layer struc- 
ture was also reported by both Eremthcenko et al.— and Kafer et al.— Furthermore, molecular 
rearrangements involving such topological phase transformations from flat to buckled pentacene 
multilayers mimicking thin film formation on Ag(l 1 1) can be expected to give the same order of 
energy differences as that of thin film and the bulk pentacene phases. One final point to be stressed 
here is that, although (i) our results may be helpful for the comparison of stability of flat and tilted 
multilayers, and (ii) the more stable multilayer configurations we found resemble the experimen- 
tally observed pentacene phases— "—'^>^ (i.e., the tilt angles are very close to that of pentacene 
bulk and thin film phases), a direct comparison of our results with the experimentally observed 
structures may not be very meaningful. This is because the multilayer and monolayer in-plane 
unit-cell dimensions which are actually different had to be chosen the same due to computational 
restrictions. 

In order to investigate the coupling of the frontier molecular orbitals of pentacene to the silver 
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Figure 5: Calculated STM image of IML Pn on Ag(lll) (a) for the occupied states around —0.9 
eV and (b)for the unoccupied states around 1.3 eV vicinity of the Fermi energy. 



substrate states, we obtained the STM pictures by using Tersoff-Hamann approximation.— The 



calculated STM images for the applied voltages of —0.9 V and 1.3 V in Figure 5 resemble to the 
HOMO and LUMO charge densities of an isolated Pn similar to Lee et a/.'s result.— Our results 
also agree well with the recent differential conductance images obtained with low-temperature 
STM experiments for seemingly similar physisorption system of pentacene/Au(lll).^^ These 



STM images in Figure 5 are consistent and are also apparent from the PDOS of IML Pn/Ag(l 11) 



presented in Figure 6 The first PDOS peak of the Pn layer about 0.5 eV below the Fermi energy 
comes from the HOMO's of the molecules. The sharpness of this peak substantiates that the fron- 
tier orbitals of the Pn molecules mixes very weakly with the 5s states of the surface Ag atoms. 
Hence, STM calculation covering an energy range of 0.9 eV below the Fermi level shows HOMO 
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charge density for Pn layer over the slab as shown in [Figure 5| where the small Ag 5s contribu 



tion was suppressed for visual convenience. Similarly, the first peak due to Pn layer about 0.7 eV 
above the Fermi energy corresponds to the LUMO's of the molecules which are also weakly mix- 
ing with the valence bands of Ag(l 11). At this point, in order to comment on the reliability of our 
method, we compare the experimental and theoretical results of Pn-Cu(lOO). Ferretti et al.,— on 
this system, have found a significant broadening of pentacene HOMO-LUMO bands and mixing 
with the substrate levels using the same computational procedure. This, in combination with the 
experimental results, were interpreted as an interaction close to chemisorption. In Pn-Ag(lll) 
case, however, it is clear from our results that the picture is more close to physisorption. 




Figure 6: (color online) Calculated PDOS for Pn/Ag(l 11) structures. The abscissa is the energy, in 
eV, relative to the Fermi level for the clean Ag(l 1 1) surface. Pentacene contributions are indicated 
by gray (red) while total DOS is in dark gray (blue). 



The bottom panel of Figure 6 shows the calculated DOS for the flat pentacene layer which is 
obtained by removing the silver substrate from 1ML/Ag(l 1 1) ( Figure 3^ -c) structure. The sharp 
peaks, having less structure, rather look like an energy level diagram due to very low overlap 
between molecular orbitals through tip to tip pentacene contacts over the layer. Metallic nature of 
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the bare Ag(l 11) surface is also presented in the succeeding panel of |Figure 6l where Fermi level is 
set as the origin of the energy axis. The DOS structure corresponding to full pentacene monolayer 



on Ag(l 1 1) is shown in the third panel from the bottom in Figure 6 DOS peaks stemming from 
pentacene show no shift in energy with respect to those of the flat Pn layer in the absence of the 
metal substrate. In addition, the broadening of each peak is localized over a small number of silver 
states indicating weak semiconductor-metal coupling. As a result of this weak interaction, the 
contact layer shows bulk-like HOMO-LUMO contributions to the total DOS around the Fermi 
energy. Therefore, STM experiments capture these frontier molecular orbitals. 

As new pentacene layers deposited on the first full monolayer the corresponding PDOS con- 
tribution starts to form localized satellite structures at around flat Pn layer peak positions. Their 
broadening is larger than the broadening in PDOS peaks obtained for 1ML/Ag(lll). This indi- 
cates that the interlayer molecular orbital overlap is relatively stronger than the coupling between 
the contact layer and the metal substrate. Moreover, these PDOS satellites match perfectly with 



the bulk pentacene DOS which is presented in the top panel of Figure 6 Therefore, energetically 
preferable thin film pentacene phase on Ag(l 1 1) up to 4ML possesses bulk- like DOS properties. 

Our DOS calculations show that pentacene has no electronic contribution at the Fermi energy. 
Evidently, highly ordered pentacene multilayers on Ag(lll), considered in this study, does not 
exhibit band transport. In addition, these multilayers occur in bulk-like phase where the overlap of 
the molecular orbitals between nearest neighbor pentacenes yield large ;r-conjugation length along 
the molecular axis. Therefore, our results suggest a hopping mechanism between the localized 
states for the carrier transport. 



Conclusion 

In conclusion, we have investigated the geometric and electronic structure of pentacene on Ag(l 11) 
surface up to 4 ML coverage at the DFT level where GGA functional perform better. At the most 
stable configuration a single isolated pentacene lies flat at 3.9 A above Ag(lll) surface on the 
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so called Bridge-60 position with a weak binding energy of —0.155 eV. For the full monolayer 
coverage, molecules align perfectly with the silver rows on the surface while the ML height extends 
slightly to 4. 1 A due to intermolecular interactions. Calculated binding energies as well as STM 
and PDOS structures indicate weak pentacene-substrate coupling. 

Pentacenes above the contact layer favor the thin film multilayer structure over the planar con- 
figuration with a slight energy difference which can be addressed to the small energy barriers 
between different phases of pentacene. Moreover, the slight misalignment of pentacene molecules 
above the first layer from the surface silver rows indicate that a bulk-like thin film phase starting 
from the second layer is adsorbed on the Ag(l 1 1) surface through a contact layer at the interface. 
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